Simultaneous cooling of axial vibrational modes in a linear ion-trap 
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In order to use a collection of trapped ions for experiments where a well denned preparation of 
vibrational states is necessary, all vibrational modes have to be cooled to ensure precise and re- 
peatable manipulation of the ions' quantum states. A method for simultaneous sideband cooling 
of all axial vibrational modes is proposed. By application of a magnetic field gradient the absorp- 
tion spectrum of each ion is modified such that sideband resonances of different vibrational modes 
coincide. The ion string is then irradiated with monochromatic electromagnetic radiation, in the 
optical or microwave regime, for sideband excitation. This cooling scheme is investigated in detailed 
numerical studies. Its application for initializing ion strings for quantum information processing is 
extensively discussed. 
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I. INTRODUCTION 

Atomic ions trapped in an electrodynamic cage allow 
for preparation and measurement of individual quantum 
systems, and represent an ideal system to investigate 
fundamental questions of quantum physics, for instance, 
related to decoherence P, , the measurement process 
0,3, or multiparticle entanglement Also, trapped 
ions satisfy all criteria necessary for quantum computing. 
Two internal states of each ion represent one elemen- 
tary quantum mechanical unit of information (a qubit). 
The quantized vibrational motion of the ions (the "bus- 
qubit" ) is used as means of communication between indi- 
vidual qubits to implement conditional quantum dynam- 
ics with two or more qubits 6]. In recent experiments 
quantum logic operations with two trapped ions were re- 
alized and teleportation of an atomic state has been 
demonstrated @. 

These implementations of quantum information pro- 
cessing (QIP) with trapped ions require that the ion 
stringis cooled to low vibrational collective excitations 
[fl I7II9I Hoj|. In particular, this condition should be ful- 
filled by all collective vibrational modes Therefore, 
in view of the issue of scalable QIP with ion traps, it is 
important to find efficient cooling schemes that allow to 
prepare vibrationally cold ion chains. 

Cooling of the vibrational motion of two ions in a com- 
mon trap po tential has been demonstrated experimen- 
tally |l2l lia. Ilil IrH (see also ^(| for a recent review). 
This is deemed to be sufficient for a quantum information 
processor which utilizes two ions at a time for quantum 
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logic operations with additional ions stored in spatially 
separated regions If more than two ions reside in a 
common trap potential and shall be used simultaneously 
for quantum logic operations, however, the task of re- 
ducing the ions' motional thermal excitation becomes in- 
creasingly challenging with a growing number of ions and 
represents a severe obstacle on the way towards scalable 
QIP with an ion chain. Straightforward extensions of 
laser-cooling schemes for one particle to many ions, like 
sequentially applying sideband cooling |lf)j | to each one 
of the modes, becomes inefficient as the number of ions 
increases, since after having cooled the last mode, the 
first one may already be considerably affected by heat- 
ing due to photon scattering and/or due to fluctuations 
of the trap potential. Therefore, it is desirable to find a 
method that allows for simultaneous and efficient cooling 
of many vibrational modes of a chain of ions. 

In this article we propose a scheme that allows for si- 
multaneous sideband cooling of all collective modes of an 
ion chain to the ground state. This is achieved by induc- 
ing position dependent Zeeman shifts through a suitably 
designed magnetic field, thereby shifting the spectrum 
of each ion in such a way that the red-sideband tran- 
sitions of each mode may occur at the same frequency. 
Thus, by irradiating the ion string with monochromatic 
radiation all axial modes are cooled. We investigate nu- 
merically the efficiency and explore implementations of 
simultaneous sideband excitation by means of laser light, 
and alternatively, by using long-wavelength radiation in 
the radio- frequency or microwave regime 0, IT^. l2fj| . 

The remainder of this article is organized as follows: 
In section [H] the cooling scheme is outlined. Numeri- 
cal investigations of the cooling efficiency are presented 
for implementations using an optical Raman transition 
(section ITTTl) and a microwave transition (section EJ- I n 
section[3possible experimental implementations are dis- 
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cussed and the cooling scheme is studied under imperfect 
experimental conditions. The paper is concluded in sec- 
tion EH 



II. THE CONCEPT OF SIMULTANEOUS 
SIDEBAND COOLING 

A. Axial vibrational modes 

We consider TV crystallized ions each of mass m and 
charge e in a harmonic trap. The trap potential has 
cylindrical symmetry around the z— axis providing strong 
radial confinement such that the ions are aligned along 
this axis |2l|. We denote by v ri v z the radial and ax- 
ial frequencies of the resulting harmonic potential, where 
Vr v z , and by z^ ' the ions classical equilibrium po- 
sitions along the trap axis. The typical axial distance 
8z between neighboring ions scales like 5z ~ £o2iV _0 
with Co = (e 2 /(47re 77Wi)) 1/3 [2l|23- For brevity, in the 
remainder of this article, the ion at the classical equilib- 
rium position z^ is often referred to as "ion j" . 

At sufficiently low temperatures the ions vibrations 
around their respective equilibrium positions are har- 
monic and the axial motion is described by N harmonic 
oscillators according to the Hamiltonian 



N 



= 22 



(4a Q + 1/2) 
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where v a are the frequencies of the chain collective modes 
and a* a and a a the creation and annihilation operators of 
a phonon at energy hv a . We denote with Q a ,P a the 
corresponding quadratures, such that [Q a ,P Q ] = ih, and 
choose the labelling convention v\ < V2 < ■ ■ ■ < vn 
whereby v\ = v z (in this article we often refer to the 
collective vibrational mode characterized by v a as "mode 



a"). The local displacement qj 



of the ion j 



from equilibrium is related to the coordinates Q a by the 
transformation 



Qj 



(2) 



where S" are the elements of the unitary matrix S 
that transforms the dynamical matrix A, characterizing 
the ions potential, such that S~ 1 AS is diagonal. The 
frequencies, v a of the vibrational modes are given by 
■y/vZ x v\ where v a are the eigenvalues of A |25j |. The 
normal modes are excited by displacing an ion from its 
equilibrium position zj ^ by an amount qj. Thus, the 
coefficients 5" describe the strength with which a dis- 
placement qj from z^ couples to the collective mode a. 

Excitation of a vibrational mode can be achieved 
through the mechanical recoil associated with the scat- 
tering of photons by the ions. This excitation is scaled 
by the Lamb-Dicke parameter (LDP) [2(j, which for a 



single ion corresponds to ^ujh/v, where uor — hk 2 /2m 
is the recoil frequency and hk the linear momentum of 
a photon. In an ion chain we associate a Lamb-Dicke 
parameter rj a with each mode according to the equation 



Va 




(3) 



Hence, if a photon is scattered by the ion at Zj , the ion 
recoil couples to the mode a according to the relation 

V" = Sf Va • (4) 

In the remainder of this article we will assume that the 
ions are in the Lamb-Dicke regime, corresponding to the 

fulfillment of condition \J (aaa a )rj a < 1. In this regime 
the scattering of a photon does not couple to the vibra- 
tional excitations at leading order in this small parame- 
ter, while changes of one vibrational quantum Kv a occur 
with probability that scales as |?y a | 2 . Changes by more 
than one vibrational quantum are of higher order and are 
neglected here. 

B. Sideband cooling of an ion chain 

In this section we consider a schematic description of 
sideband cooling of an ion chain, in order to introduce the 
concepts relevant for the following discussion. We denote 
by |0) and |1) the internal states of the ion transition at 
frequency loq, in absence of external fields, and linewidth 
7. A spatially inhomogeneous magnetic field is applied 
that shifts the transition frequency of each ion individu- 
ally such that for the ion at position z^ ' the value ujj is 
assumed. Each ion transition couples to radiation at fre- 
quency u>l, which drives it well below saturation. In this 
limit, the contributions of scattering from each ion to the 
excitation of the modes add up incoherently [53, ■ 

For this system, the equations describing the dynamics 
of laser sideband cooling of an ion chain can be reduced 
to rate equations of the form 
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1) A a _P a {n {a ^ + I) - A%P a (n {a ^) 
A a _P a {n^)- A%P a (n^ -1)1 (5) 



where P a irS°^) is the average occupation of the vi- 
brational number state \rS a ^) of the mode a, and A" 
(A^_ ) characterizes the rate at which the mode is heated 
(cooled). Equation (J5J is valid in the Lamb-Dicke regime, 
i.e. when the LDP is sufficiently small to allow for a per- 
turbative expansion in this parameter. Denoting by ttj 
the Rabi frequency, the heating and cooling rate takes 
the form El 



A% 



N 



n 2 

27 
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(6) 
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where the detuning Sj = uji 



The coefficient 



emerges from the integral over the angles of photon emis- 
sion, according to the pattern of emission of the given 
transition j2||. For At > A^_ a steady state exists, it is 
approached at the rate 

-,(<*) 



1 cool — A - 



(7) 



cool — + V / 

and the average number of phonons of mode a at steady 
state is given by the expression 



A% 



A% 



(8) 



Sideband cooling reaches (rv- a >) <C f through At S> At- 
This condition is obtained by selectively addressing the 
motional resonance at u>q — u a . This is accomplished for 
a single collective mode when 7 -C v a and 5 a = v a . 

In this work, we show how the application of a suit- 
able magnetic field allows for simultaneous sideband cool- 
ing of all modes. In particular, the field induces space- 
dependent frequency shifts that suitably shape the ex- 
citation spectrum of the ions. Simultaneous cooling is 
then achieved when for each mode a there is one ion j 
with the matching resonance frequency, that is, such that 



dj = UJ L 



This procedure is outlined in detail 



in the following subsection. 



C. Shaping the spectrum of an N ion chain 
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FIG. 1: Illustration of the axial motional spectrum of a chain 
of 10 ions in the presence of a spatially inhomogeneous mag- 
netic field. The vertical lines indicate the axial position of 
each ion in units of (q. The corresponding horizontal lines 
indicate the frequencies of the spectral lines as measured at 
that particular ion. The thick horizontal lines indicate the 
ions resonance frequencies uij — uji (in units of the secular ax- 
ial frequency v\) relative to the resonance frequency uji of the 
ion at z\ — — 2.87£o. The remaining horizontal lines show the 
frequencies of red sideband resonances for each ion at frequen- 
cies uij — v a (j, a = 1, . . . , 10). The magnetic field is designed 
such that — v\ — U12 — ^2 = • • • wjv — vn ■ These resonances 
are highlighted by medium thick lines. 



Assume the ion transition |0) — > |1) and that a mag- 
netic field-whose magnitude varies as a function of z-is 
applied to the linear ion trap, Zeeman shifting this res- 
onance. As a result, the ions resonance frequencies Wj 
are no longer degenerate. The field gradient is designed 
such that all ions share a common motion-induced reso- 
nance. This resonance corresponds to one of the transi- 
tions |0, n^) — > \1, 7i/ Q ) — 1), namely to the red sideband 
of the modes a. The resonance frequency of each ion is 
shifted such that the red sidebands of all modes can be 
resonantly and simultaneously driven by monochromatic 
radiation at frequency to = ijQ\ — v\ = . . . = uj^—vn- Ionic 
resonances and the associated red sideband resonances 
optimally shifted for simultaneous cooling are illustrated 
in Fig. n for the case of 10 ions. 

Sideband excitation can be accomplished by either 
laser light or microwave radiation according to the 
scheme discussed in ^8|. With appropriate recycling 
schemes this leads to sideband cooling on all N modes 
simultaneously. A discussion on how a suitable field gra- 
dient shifting the ionic resonances in the desired fashion 
can be generated is deferred to section Ivl 

D. Theoretical model 

As an example, we discuss simultaneous sideband cool- 
ing of the collective axial modes of a chain composed of 




FIG. 2: Schematic of the ions internal energy levels on which 
cooling is implemented. Indicated are the relevant Rabi fre- 
quencies (symbols Q), spontaneous decay rates (r), and de- 
tunings (A). The corresponding equations for the dynamics 
are discussed in detail in the appendix 



1 Yb + ions with mass m = 171 a.m.u.. The ions are 
crystallized along the axis of a linear trap characterized 
by v\ — 1 x 27rMHz. A magnetic field B(z) along the axis 
is applied that Zeeman-shifts the energy of the internal 
states. The value of the field along z is such that it shifts 
the red-sidebands of all modes into resonance along the 
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chain, while at the same time its gradient is sufficiently 
weak to negligibly affect the frequencies of the normal 
modes [lflj . 

The selective drive of the motional sidebands can 
be implemented on a magnetic dipole transition in 
m Yb + close to ujq = 12.6 x 2n GHz between the hy- 
perfine states |0) = \S 1/2 ,F = 0) and |1) = \S 1/2 ,F = 
l,mp = 1). The magnetic field gradient lifts the degen- 
eracy between the resonances of individual ions, and the 
transition frequency u>j of ion j is proportional to B(zj) 
in the weak field limit [AbB/Hujo <C 1, where [1b is the 
Bohr magneton. For strong magnetic fields the variation 
of ujj with B is obtained from the Breit-Rabi formula Q . 

We investigate two cases, corresponding to two dif- 
ferent implementations of the excitation of the sideband 
transition between states |0) and In the first case, 
discussed in section ITTT1 the sideband transition is driven 
by two lasers with appropriate detuning, namely a Ra- 
man transition is implemented with intermediate state 
|2) = \P\/2)- In the second case, presented in section llVl 
microwave radiation drives the magnetic dipole. 

Since spontaneous decay from state |1) back to |0) is 
negligible on this hyperfine transition, laser light is used 
to optically pump the ion into the |0) state via excitation 
of the |1) — > 1 2) electric dipole transition. This laser light 
is close to 369nm and serves at the same time for state se- 
lective detection by collecting resonance fluorescence on 
this transition, and for initial Doppler cooling of the ions. 
The state |2) decays with rates I^i = 11 x 27rMHz and 
= 5.5 x 27rMHz into the states |1) and |0), respec- 
tively 29]. The considered level scheme is illustrated in 
Fig. [21 and the corresponding model is described in the 
appendix. 

We evaluate the efficiency of the cooling procedure 
by neglecting the coupling between different vibrational 
modes by photon scattering, which is reasonable when 
the system is in the Lamb-Dicke regime. In this case, 
the dynamics reduce to solving the equations for each 
mode a independently, and the contributions from each 
ion to the dynamics of the mode are summed up inco- 
herently [2]j, as outlined in Sec. Ill Bl The steady state 
and cooling rates for each mode are evaluated using the 
method discussed in and extended to a chain of N 
ions. The extension of this method to a chain of ions is 
presented in the appendix. The numerical calculations 
were carried out for this scheme and chains of N ions 
with 1 < N < 10 and for some values N > 10. Since the 
qualitative conclusions drawn from these calculations did 
not depend on TV, we therefore restrict the discussion in 
sections imiUVl and El to the case N = 10. 



III. RAMAN SIDEBAND COOLING OF AN ION 
CHAIN 

We consider sideband cooling of an ion chain when the 
red sideband transition is driven by a pair of counter- 
propagating lasers, which couple resonantly the levels 



|0) and |1). The two counter-propagating light fields 
couple with frequency ojri, lor2 to the optical dipole 
transitions |0) — > |2) and |1) — > |2), respectively. The 
two lasers are far detuned from the resonance with level 
1 2) such that spontaneous Raman transitions are negligi- 
ble compared to the stimulated process. We denote by 
A01 = [{uri —urz) — the Raman detuning, such that 
A01 = corresponds to driving resonantly the transition 
|0) — > |1) at the first ion in the chain, and by f2oi the 
Rabi frequency describing the effective coupling between 
the two states. A third light field with Rabi frequency 
SI12 is tuned close to the resonance |1) — > |2) and serves 
as repumper into state |0) (compare Fig.|2). The frequen- 
cies uJm are close to the 171 Yb + resonance at 369nm, and 
the trap frequency is v = 1 x 27rMHz. Hence, from Eq. @ 
the Lamb-Dicke parameter takes the value 771 ss 0.0926. 



A. Sequential cooling 

In absence of external field gradients shifting inhomo- 
geneously the ions transition frequencies (namely, when 
lui = ... = lun = ujq), cooling of an ion chain could 
be achieved by applying sideband cooling to each mode 
sequentially. In each step of the sequence all ions are 
illuminated simultaneously by laser light with detuning 
A01 = — v a , thereby achieving sideband cooling of a par- 
ticular mode a. Since all ions are illuminated, they all 
contribute to the cooling of mode a. 

In Fig.|3K) the steady state vibrational number of each 
mode at the end of the cooling dynamics is displayed as 
a function of the relative detuning A i. Each mode v a 
reaches its minimal excitation at values of the detuning 
A01 = — v a . Therefore, in order to cool all modes close 
to their ground state, the detuning of the laser light has 
to be sequentially set to the optimal value for each mode 
a. 

The cooling rates of mode a at A01 = — v a , as de- 
fined in Eq. J7J, arc displayed in Fig. [3)). They arc differ- 
ent for each mode and vary between 1kHz and 100kHz for 
the parameters chosen here. Even though these cooling 
rates would, in principle, allow for cooling sequentially all 
modes in a reasonably short time, this scheme may not 
be effective, since while a particular mode a is cooled all 
other modes are heated (i) by photon recoil, and, (ii) by 
coupling to the environment. As external source of heat- 
ing we consider here the coupling of the ions charges to 
the fluctuating patch fields at the electrodes [32|. The 
effects of these processes on the efficiency of cooling are 
discussed in what follows. 

The consequences of heating due to photon scattering 
are visible in Fig.[3Ji). Here, one can see that while cool- 
ing one mode, others can be simultaneously heated, such 
that their average phonon number at steady state is very 
large. These dynamics are due to the form of the reso- 
nances in a three- level configuration [3(], |3^| . In general, 
however, the time scale of heating processes due to pho- 
ton scattering is considerably longer than the time scale 
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FIG. 3: Raman sideband cooling of a chain of 10 ions with- 
out magnetic field gradient. The parameters are Qi 2 = 
100 x 27rkHz, f2 i = 30 x 2nkRz, A12 = -10 X 2ttMHz. a) 
Steady state mean vibrational number (n^"') as a function 
of A01 in units of v\. b) Cooling rate F^j of mode a for 
A i = — v a . c) Steady state populations for A i = —Vi, 
corresponding to sideband cooling of mode 1. The bars in- 
dicating the vibrational excitations of modes 1/3 and V4 have 
been truncated. 



at which a certain mode is optimally sideband cooled, 
since the transitions leading to heating are out of reso- 
nance. In the case discussed in Fig. [31 for instance, the 
heating rates of these modes at A i = —v\ are orders 
of magnitude smaller than the cooling rate of mode 1, 
and their dynamics can be thus neglected while mode 
1 is sideband cooled. Similar dynamics are found for 
Aqi = — v a . Thus, in general one may neglect photon 



scattering as source of unwanted heating of modes that 
are not being efficiently cooled. 

Nevertheless, heating by fluctuating electric fields oc- 
curs with appreciable rates ranging between 5s _1 and 
10 4 s^ 1 [33. The heating rate is different for each mode 
and was observed to be considerably larger for the COM 
mode (here denoted as mode 1) than for modes that 
involve differential relative displacements of individual 
ions. Obviously, cooling can only be achieved, if the cool- 
ing rate soq °f eacn mode exceeds in magnitude the 
corresponding trap heating rate denoted by r h " 



heat ' 



p(«) ^ p(«) w _ 1 

1 heat 1 cool,seq vu! — x > 
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(9) 



In addition, one must consider that after a particular 
mode a has been cooled, it might heat up again while all 
other modes, /3 ^ a are being cooled. This imposes a sec- 
ond condition on the cooling rate. In order to quantify 
this second condition, we first evaluate the time, T r 



(a) 
cool 



it takes to cool one particular mode a from an initial 
thermal distribution, obtained by means of Doppler cool- 
ing and characterized by the average occupation number 



to a final distribution characterized by 



This time can be estimated to be |2£ 



CODl (n(*>) 



/(r (a) 

i V cool.scq 



heat/ 



(10) 



where (n^) / < 1 at steady state was assumed. 



From this relation one obtains the total time, T E 



(a) 



SCCj 



needed to cool all modes except mode a, or, in other 
words the time during which mode a is not cooled and 
could get heated. This time, Tshq needed to sideband 
cool all modes with (3 ^ a is 



71(a) = T (0) 

seq / j cool 



i/(r 



C/3) _ p(/3) 
cool,scq heat 



)• 



(11) 

If mode a is to stay cold during this time, the heating rate 
affecting it must be small enough. Hence the condition 
for efficient sequential sideband cooling is derived, 



p( Q ) x 71(a) ^ 1 

1 heat X 1 seq ^ 1 
.(a) 



(12) 



namely, during time T seq , necessary for cooling the 
modes (3 ^ a, the heating of mode a has to be negligible. 
Clearly, this condition is stronger than the one derived 
in relation © , and its fulfillment becomes critical as the 
number of vibrational modes (ions) is increased. 

A rough estimate of the time T^q to be inserted in (fT2|) 
can be obtained from eq. (|llfl under the assumption that 
all modes start out with the same mean excitation (n)j 
(usually determined by initial Doppler cooling) and are 
cooled to the same final excitation (n) / . Using condi- 
tion ©, one obtains 



71(a) 

seq 



In 



(n)i 



cool.scq 



(13) 
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Substituting this expression into (|12f> gives 

itL « fin £ ^ r<, scq (14) 

which places a stronger restriction than © on the trap 
heating rate that can be tolerated, if sequential cooling is 
to work. This relation has to hold true for a = 1, . . . , N. 
Expression l|14|) will be used for a comparison with si- 
multaneous sideband cooling (see Sec. IIII C|l . 

B. Simultaneous cooling 

We consider now the case, when a magnetic field gra- 
dient is applied to the ion chain, such that the situation 
shown in Fig. ^ is realized. The axial modes of the chain 
can then be simultaneously cooled. 

Figure0]displays the steady state mean vibrational ex- 
citations that are obtained when the effective Rabi fre- 
quency for the Raman coupling f2oi = 5 x 27rkHz, the 
Rabi frequency of the repumper f2i2 = 100 x 27rkHz, 
and the detuning A 12 = -10 x 27rMHz. Fig. dis- 
plays the mean vibrational quantum number (n^) of 
the COM mode as a function of the detuning Aqi. Here 
three minima are visible. The leftmost minimum occurs 
at Aoi = —v\ and corresponds to resonance with the 
red sideband of the COM in the spectrum of the first 
ion. The minimum in the middle stems from the reso- 
nant drive of the red COM-sideband in the spectrum of 
the second ion while the one on the right is caused by the 
spectrum of the third ion in the chain. The location of 
these resonances correspond to the ones shown in Fig. ^ 
Heating of the COM mode occurs if the blue sideband of 
the COM mode is driven resonantly. In Fig.^Ji) the heat- 
ing at the blue sideband of the first ion, i.e. Aoi = v\ , is 
visible. 

Figure^)) displays (n^), (ri (2) ), (n 1 -^), and (n (4) ) as 
a function of Aoi. These mean excitations have been cal- 
culated using the same parameters as in Fig. 0]}) . The 
minima visible in this figure can be identified with the 
corresponding resonances in the spectra of the ions by 
comparison with Fig. ^ A common minimum occurs at 
Aqi = — v\ where all four vibrational modes are simulta- 
neously cooled to low excitation numbers. 

The mean vibrational quantum number (nA a )) of all 
ten axial modes is displayed in Fig. |3J;) as a function 
of the detuning Aoi in the neighborhood of the value 
Aoi = — v\. At this value of Aoi the mean excitation 
(n^O) reaches its minimum for all modes. The mode 
at frequency z^io displays a relatively narrow minimum 
and its mean vibrational number, although very small 
at exact resonance, is orders of magnitude larger than 
the ones of the other modes. In fact, ion 10 participates 
only little in the vibrational motion of mode 10. This 
is described by the small matrix element S\q — 0.0018 
that scales the corresponding Lamb-Dicke parameter as 
shown in Eq. (gj. 
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FIG. 4: Raman sideband cooling a chain of 10 ions in the pres- 
ence of a spatially inhomogeneous magnetic field such that the 
condition u>i — V\ = ... = ujn — vn is fulfilled (see Fig. 0. 
The parameters are Q.12 = 100 x 27rkHz, J7oi = 5 x 27rkHz, 
A12 = — 10 x 27rMHz. The steady state mean vibrational 
excitation (n' a '} is displayed, a) (n' 1 ') (COM mode) as a 
function of Aoi in units of vi. b) (n^ a ') with a = 1,2,3,4 
as function of Aoi. c) {ro' a ') for a = 1, . . . , 10 as a func- 
tion of Aoi. d) (n' a ') for a = 1, . . . , 10 when the chain is 
simultaneously cooled at the detuning Aoi = — V\. 
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Fig. 0Ji) displays the steady state temperature of each 
mode when the detuning of the Raman beams is set close 
to —v\. At this detuning the average excitation reaches 
its minimum for each mode, which is (n^) < 10 -3 . 



Cooling rates for simultaneous Raman cooling 



5 6 
Mode 



FIG. 5: Cooling rates for Raman sideband cooling a chain of 
10 ions in the presence of a spatially inhomogeneous magnetic 
field such that the condition u/i — v\ — ... = cum — vn is 
fulfilled. The same parameters have been used here as for 
generating Fig. [3 (fi 12 = 100 x 27rkHz, Q i = 30 x 27rkHz, 
A12 = — 10 x 2-7rMHz). The black bars indicate the rate for 
each mode when simultaneously sideband cooling all modes. 
The grey bars give the cooling rates that are achieved, if (still 
in the presence of the magnetic field) the Raman detuning is 
set to that value where the maximum cooling rate for each 
individual mode is obtained. 



We now turn to the cooling rates that are achieved 
when simultaneously cooling all modes, that is, the field 
gradient leading to the spectrum in Fig.l is applied, and 
Aoi = —V\. These rates are indicated by black bars 
in Fig. [5] and have been evaluated with the same pa- 
rameters used for the simulation of sequential cooling in 
Fig. El that is, Q i = 30 x 27rkHz, Q 12 = 100 x 27rkHz, 
and A12 = -10 x 2ttMHz. 

Figure [3] shows that the COM mode characterized by 
V\ and the mode characterized by v 2 are cooled most 
efficiently, while modes 9 and 10 display much smaller 
cooling rates. In particular, » 10 _1 s _1 which will 

make cooling of this mode very slow at Aoi = — v\ |4l|. 

The origin of this behaviour can be understood as fol- 
lows. The scheme of simultaneous sideband cooling re- 
quires that each ion of the chain is employed to cool one of 
the axial modes. This is achieved by applying a suitable 
magnetic field gradient. The simplest experimental im- 
plementation uses a monotonically increasing magnetic 
field, such that the first ion of the chain is used for cool- 
ing mode 1, the second ion mode 2, etc. (compare Fig. 
1). However, the coupling of a certain ion displacement 
to a certain mode can be very small. It occurs, for in- 
stance, that the largest axial excitations couple weakly to 



the ions at the edges of the ion string, instead they are 
mainly characterized by oscillations of ions in the center 
of the chain [34[ . A manifestation of this behavior is the 
small value of the matrix element Si®. Thus, mode 10 is 
not efficiently cooled by illuminating ion 10, rather it is 
optimally cooled by addressing an ion closer to the cen- 
ter of the chain. This is evident by inspection of the grey 
bars in Fig. [3] that indicate the optimal cooling rate for 
each individual mode, obtained by employing that par- 
ticular ion in the chain which has the largest coupling to 
the mode to be cooled. In presence of the magnetic field 
gradient, this optimal cooling rate is achieved, if the de- 
tuning, Aoi is set such that the appropriate red-sideband 
transition of this particular ion is driven resonantly. In 
this way, the cooling rate is maximal for that particu- 
lar mode (i.e., each grey bar corresponds to a different 
detuning) . 

For the case of mode 10 the difference between the 
cooling rates depending on which ion is addressed is par- 
ticularly striking: the cooling rate of mode 10 at detun- 
ing Aoi = —v\ (black bar) is very small, as noted above, 
whereas the optimal rate = 2.63 x 10 3 s _1 (given the 
parameters used here) is achieved at A i = —3.28 x v\. 
At this detuning ion 6 is used for cooling mode 10, as can 
be seen from Fig. 

Efficient cooling of all modes, given the field gradi- 
ent that gives rise to the spectrum illustrated in Fig. 1, 
can be obtained by combining sequential and simultane- 
ous cooling as follows: First the ion string is illuminated 
with radiation such that Aoi = —3.28 x v\ which will 
efficiently cool mode 10 and will have little effect on all 
other modes. Then, we apply radiation with Aoi = —v\ 
(cooling rates indicated by black bars in Fig. [SJ. This 
will simultaneously cool all other modes. This is a spe- 
cific recipe to cool 10 ions. For an arbitrary number of 
ions, first one cools the modes which cannot be efficiently 
simultaneously cooled with all other modes, and then, in 
a second step, as many modes as possible are cooled si- 
multaneously (in our case modes 1 through 9). 

A necessary condition for the efficiency of this scheme 

is 



p(«) ^ p( Q ) 

heat cool.sim 



V a = 1, 



10 . 



(15) 



where Y 



(a) 

COOl,! 



is the cooling rate of mode a when si- 



multaneously cooling it with the other modes. Condi- 
tion (|15fl is the analogue to 10 which we derived for se- 
quential cooling. Moreover, since in the scheme proposed 
here modes 1 through 9 are simultaneously cooled, a con- 
dition analogous to the one in Eq. (|12|) has to be fulfilled 
for mode 10 only: After mode 10 has been cooled, this 
mode may not heat up again appreciably during the time 
for which modes 1 through 9 are simultaneously cooled. 
This condition is expressed as 



.(10) 
heat 



(16) 



where is the time during which modes 1 through 

9 are cooled simultaneously to the desired final values, 
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that is, the time during which mode 10 is not cooled and 
could heat up. Here, it is given by 



T s \m = max /3 =i...9 



In- 



/(r 



(/?) 
cool 



r (/3) \ 
A heat-' 



(17) 



This relation is analogous to relation l|ll|) for sequential 
cooling. 



The times T^q (eq. (DJ) and T sim ( ec l- C3) S ive an 
upper limit for the tolerable trap heating rate of mode 
10 for the cases of sequential and simultaneous cooling, 
respectively. It turns out that the relevant time scales 
Tseq^ and are of the same order of magnitude: Us- 

ing eq. El to compute for the parameters used 

in Fig. gives « ln((n)i/(n)/) x 2.4 x 10~ 3 s, 

while Tscq' 1 ~ lix((n)i/{n) t) x 1.3 x 10~ 3 s is obtained 
from eq. (|11|1 for the same set of parameters j^. Thus, 
the cooling scheme proposed here does not increase the 
admissible trap heating rate of mode 10. Since this rate 
is expected to be low anyway, this will not be a restric- 
tion prohibiting the successful cooling of a long ion chain 
using either sequential or simultaneous cooling. 



C. Comparison of sequential and simultaneous 
cooling 

So far we have stated the general conditions that have 
to be met in order to efficiently cool an ion chain with 
sequential cooling and with a scheme combining simul- 
taneous and sequential sideband cooling. In this section 
we discuss their efficiencies. 

We note that when using sequential sideband cooling, 
one may utilize all ions in the chain in order to cool one 
mode, where the cooling rates of each ion add up inco- 
herently. In the case of simultaneous sideband cooling, 
on the other hand, only one ion is employed in order to 
cool a particular mode. Assuming that the coupling of 
this ion to the mode is sufficiently large to allow for ef- 
ficient cooling, the following expression for the average 
simultaneous cooling rate is deduced 



r ~lr 

1 cool, sim ~ cool.scq 



(18) 



where r coo i iSeq is the cooling rate for sequential cooling 
averaged over all modes. Thus, for simultaneous cooling 
relation l|15|) yields 



L heat j\f cool > ; 



r< 



(19) 



From (|18fl . since the total cooling rate is r <scq = 

J- cool.scq 

/N, one obtains that the efficiencies of simulta- 
neous and sequential cooling are comparable. However, 
it should be remarked that this estimate corresponds to 
the worst case for simultaneous cooling. In fact, esti- 
mate 1)181 is correct for long- wavelength vibrational exci- 
tations, which correspond to low-frequency axial modes, 



where practically all ions of the chain participate in the 
mode oscillation. Short-wavelength excitations, on the 
other hand, are characterized by large displacements of 
the central ions, while the ions at the edges practically 
do not move j34j. Due to this property, for these modes 
one may find a magnetic field configuration such that 

p(a) ^ p(a) 

cool, sim cool,scq' 

On the basis of these qualitative considerations, one 
may in general state that simultaneous cooling of the 
chain is at least as efficient as sequential cooling. For 
the configuration discussed in this paper the two meth- 
ods are comparable. Differently from sequential cooling, 
however, the efficiency of simultaneous cooling can be 
substantially improved by choosing a suitable magnetic 
field configuration that maximizes coupling of each mode 
to one ion of the chain. 



IV. SIDEBAND COOLING USING 
MICROWAVE FIELDS 

A. Effective Lamb-Dicke parameter 

We investigate now sideband cooling of the ion chain's 
collective motion using microwave radiation for driving 
the sideband transition. It should be remarked that 
in this frequency range sideband excitation cannot be 
achieved by means of photon recoil which for long wave- 
lengths is negligible. This is evident from eq. fusing a 
typical trap frequency, v a of the order 27r x 1MHz. Never- 
theless, in 0| it was shown that with the application of 
a magnetic field gradient an additional mechanical effect 
can be produced accompanying the absorption/emission 
of a photon. This is achieved by realizing different me- 
chanical potentials for the states |1) and |0) P. ITsl IT^|. 
Hence, by changing an ion's internal state by stimulated 
absorption or emission of a microwave photon the ion 
experiences a mechanical force. 

An effective LDP can be associated with this force, 
that is defined as 1131: 



'ha ° 



IE A 



where 



£< 3 ~ S 3 



a d z ujj Az a 



(20) 



(21) 



Here d z ojj is the spatial derivative of the resonance fre- 
quency of the ion at z]p with respect to zj, and Az a — 

y/h/ (2mv a ). The reader is referred to the appendix for 
the theoretical description of the ion chain dynamics in 
presence of a spatially- varying magnetic field. 

The term 77" appearing in Eq. (|20|l is the LDP due to 
photon recoil and defined in Eq. while the term e" 
is the LPD arising from the mechanical effect induced by 
the magnetic field gradient. Their ratio is given by 



kAz n 



2tt Az n Kj 



(22) 
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where A is the wavelength of the considered transition 
and = Az a d z LUj/i> a is the rescaled frequency gradient 

Mm 

For a transition in the microwave frequency range, like 
the transition between the states |0) and |1) in 171 Yb + , 
using typical values, like k" w 10~ 3 , A w 10~ 2 m, Az a w 
10 _8 m, one finds that the Lamb-Dicke parameter due to 
the recoil of a microwave photon is at least two orders of 
magnitude smaller than the one due to the mechanical 
effect induced by the magnetic field, and thus ?y°f ss e" 
in the microwave region. 

B. Steady state population and cooling rates 

We study sideband cooling using microwave radiation 
using the model outlined in the appendix. For the Rabi 
frequencies the same parameters as in section IIII Bl are 
employed (i.e., f2 i = 5 x 27rkHz, f2i 2 = 100 x 2-7rkHz, 
A 12 = -10 x 2ttMHz). The effective LDP is determined 
by the magnitude of the magnetic field gradient that is 
used to superimpose the motional sidebands of all vi- 
brational modes, and is given by rf-^ ~ e" = Sfn" ~ 
Sf x2x 10~ 3 . 

It must be remarked that these parameters are out- 
side the range of validity for the application of the nu- 
merical method we use. In fact, by applying it one ne- 
glects the fourth and sixth orders in the LDP expansion 
of the optical transition of the repumping cycle, which are 
of the same order of magnitude as the microwave side- 
band excitation. Nevertheless, given the complexity of 
the problem, characterized by a large number of degrees 
of freedom, we have chosen to use this simpler method 
in order to get an indicative estimate of the efficiency as 
compared to the case in which the sidebands are driven 
by optical radiation. Therefore, the results we obtain in 
this section are indicative. In fact, neglecting higher or- 
ders in the LDP of the optical transition corresponds to 
underestimate heating effects due to diffusion. 

Fig. EJ0 shows the steady state mean vibrational quan- 
tum number (n^) f of the 10 axial vibrational modes as a 
function of the detuning Aoi of the microwave radiation, 
where Aoi = when the microwave field is at frequency 
u>i. Figure^) displays the final excitation number of all 
vibrational modes as a function of Aoi around the value 
Aoi = Figure shows the steady state mean 

excitation of all 10 modes at Amw = —V\. 

All vibrational modes are cooled close to their ground 
state. However, the average occupation number of the 
highest vibrational frequency v\q is orders of magnitude 
larger than the ones of the other modes. The origin of 
this behavior is the small value of the coefficient S[q — 
0.0018 as is discussed in Sec. 1111 Hi A comparison of Fig. 

with the results obtained by simultaneously cooling 
using an optical Raman processes (Fig. 0Ji) shows that 
the mean vibrational numbers that are achieved here are 
considerably larger. 

In Fig. EJi) the cooling rates are displayed that are 



a) 




-7 -6 -5 -4 -3 -2 -1 
A„ w (untt5QfvJ 




123456789 10 
Mode 



dj 




1 23456789 10 



FIG. 6: Microwave sideband cooling a chain of 10 ions in 
the presence of a spatially inhomogeneous magnetic field such 
that the condition ui\ — v\ = . . . = ujm ~ vn is fulfilled, that is, 
simultaneous sideband cooling using microwave radiation for 
sideband excitation. Same parameters as for Fig. 0] (f2i2 = 
100 x 27rkHz, f2 i = 5x27rkHz, A i2 = -10x2ttMHz). a) Mean 
vibrational numbers (n'™')/ at steady state as a function of 
the detuning Aoi • b) Same as a) for the values around Aoi = 
—v\. c) {rv' c '')f at Aoi = —vi, i.e. when all ten vibrational 
modes are simultaneously cooled, d) Cooling rates for Aoi = 
—v\. 
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obtained when driving the sideband transition with mi- 
crowave radiation. These rates are much smaller than the 
ones shown in Fig. that are obtained using an optical 
Raman process. 

From this comparison it is evident that simultaneous 
cooling is more effective by using an optical Raman pro- 
cess than microwave radiation. In particular, the cool- 
ing rates obtained with simultaneous microwave sideband 
cooling can be comparable to the heating rates in some 
experimental situations, resulting in inefficient cooling. 

This striking difference in the efficiency can already be 
found, if comparing the two methods when they are ap- 
plied to cooling a single ion. Its origin lies in the fact 
that in the optical case the LDP accounting for the pho- 
ton recoil, namely 77™ in Eq. I|20|) . is considerably larger 
than the LDP, e" caused by the magnetic field gradi- 
ent used here to superimpose the motional sidebands. 
This can be verified by using an optical wavelength in 
Eq. (12211 which gives r]°^ ~ iff . On the other hand, driv- 
ing the |0) — |1) transition directly by microwaves results 
in jf^ w s" <C iff ■ For the model system considered 
here, the optical LDP is at least one order of magnitude 
larger than the microwave LDP for simultaneous cooling. 

In the microwave sideband cooling scheme optical tran- 
sitions are used for repumping, thus leading to an en- 
hanced diffusion rate during the dynamics and thus to 
lower efficiencies. The fundamental features of these dy- 
namics can be illustrated by a rate equation of the form 
Eq. |J5J , describing sideband cooling of a single ion with 
microwave radiation, where the rates are 
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(23) 
,(24) 



ffopt accounts for the recoil due to the spontaneous emis- 
sion when the ion is optically pumped to the state |0) |27| . 
and e is the LDP for the microwave transition due to 
the field gradient. The latter multiplies the terms where 
a sideband transition occurs by microwave excitation, 
whereas ?7 pt multiplies the terms where sideband exci- 
tation occurs by means of spontaneous emission, which 
thus describe the diffusion during the cooling process. 
Efficient ground state cooling is achieved when the rate 
of cooling is much larger than the rate of heating, which 
corresponds to the condition A-/A+ ^> 1. For v\ ^> 7^, 
whereby 7^ is the linewidth of the |0) — > |1) transition, 
this ratio scales as 



A. 
~A, 



Av 2 



Iffoptl 2 7 2 



(25) 



This result differs for the ratio obtained in the all-optical 
case, where A-/A+ ~ Av 2 /j 2 . For typical parameters, 
corresponding to a magnetic field gradient that super- 
poses all sidebands in an ion trap, |e| <C |f?opt|- Hence, 
for a given value of the ratio j/v the cooling efficiency 
in the optical case is considerably larger than in the mi- 
crowave case. 



Note that the parameter e can be made larger by in- 
creasing the magnitude of the magnetic field gradient, 
which in Eq. H22[l corresponds to increasing k". However, 
if the modes of an ion chain are to be cooled simultane- 
ously, the choice of the magnetic field gradient is fixed by 
the distance between neighboring ions, and the efficiency 
is thus limited by this requirement. 

If an individual ion (or a neutral atom confined, for ex- 
ample, in an optical dipole trap) is to be sideband cooled, 
or sequential cooling is applied to a chain of ions, then 
the above mentioned restriction on the magnitude of e is 
not present and microwave sideband cooling can be as ef- 
ficient as Raman cooling. For this case, method [3(j may 
be implemented, provided the Lamb-Dicke regime applies 
and ?7opt and e are of the same order of magnitude. 



V. EXPERIMENTAL CONSIDERATIONS 

In this section we discuss how the magnetic field gradi- 
ent for simultaneous sideband cooling can be generated 
and how cooling is affected, if the red motional sidebands 
of different vibrational modes are not perfectly super- 
posed. In order to demonstrate the feasibility of the pro- 
posed scheme it is sufficient to restrict the discussion to 
very simple arrangements of magnetic field generating 
coils. 

The use of aposition dependent ac-Stark shift has been 
proposed in |36( to modify the spectrum of a linear ion 
chain. This may be another way of appropriately shifting 
the sideband resonances but will not be considered here. 



A. Required magnetic field gradient 

If the vibrational resonances and the ions were equally 
spaced in frequency and position space, respectively, then 
a constant field gradient, appropriately chosen, could 
make all N modes overlap and let them be cooled at 
the same time. Since v a — v a -\ decreases monotonically 
with growing a and the ions' mutual distances vary with 
Zj, the magnetic field gradient has to be adjusted along 
the z— axis. The field gradient needed to shift the ions' 
resonances by the desired amount is obtained from 



dB 

dz 



(z j +z j - 1 )/2 



B{zj) - B( Zj -{) 

Zj - Zj-t 

G - G-i mb 



•,(86) 



where Q = Zj/Co is the scaled equilibrium position of 
ion j, and Vj is the square root of the j— th eigenvalue 
of the dynamical matrix. Eq. 1261 describes the situa- 
tion for moderate magnetic fields (the Zeeman energy 
is much smaller than the hyperfine splitting), such that 
d z ujj = 1/2 g.j^Bd z B with gj w g s = 2 (state |0) does 
not depend on B). As an example, we consider again a 
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string of N = 10 m Yb + ions in a trap characterized by 
v z = lx 27rMHz (thus, Co = 2.7/im). 
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FIG. 7: Required magnetic field gradient to superpose the 
motional red sidebands of ten 171 Yb + ions (markers) in a 
trap characterized by a COM frequency ci=lx 27rMHz and 
calculated field gradient (solid line) produced by three single 
wire windings (see text). 



The markers in Fig. [7] indicate the values of the re- 
quired field gradient according to eq. 1261 whereas the 
solid line shows the gradient generated by 3 single wind- 
ings of diameter 100/im, located at z = —100,50, and 
100 \im w 36Co, respectively (the trap center is chosen as 
the origin of the coordinate system) ■ Running the 
currents -5.33A, -6.46A, and 4. 29A, respectively, through 
these coils produces the desired field gradient at the lo- 
cation of the ions. Micro electromagnets with dimen- 
sions of a few tens of micrometers and smaller are now 
routinely used in experiments where neutral atoms are 
trapped and manipulated |38|. Current densities up to 
10 8 A/cm 2 have been achieved in such experiments. A 
current density more than two orders of magnitude less 
than was achieved in atom trapping experiments would 
suffice in the above mentioned example. 

This configuration of magnetic field coils shall serve as 
an example to illustrate the feasibility of the proposed 
cooling scheme in what follows. It will be shown that 
with such few current carrying elements in this simple 
arrangement one may obtain good results when simul- 
taneously sideband cooling all axial modes. More so- 
phisticated structures for generating the magnetic field 
gradients can of course be employed, making use of more 
coils, different diameters, variable currents, or completely 
different configurations of current carrying structures. 

Ideally, all 10 sideband resonances would be super- 
imposed for optimal cooling. The resonances shown in 
Fig. [S] result from the field gradient calculated using 
the simple field generating configuration described above, 
and do not all fall on top of each other. Neverthe- 
less, Fig. [5] shows how well all 10 sideband resonances 
are grouped around u>x — v\. Vertical bars indicate the 
location relative to lo\ of the red sideband resonance, 
ujj — Vj of the j— th ion, with j = 1, . . . , 10. These 
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FIG. 8: The vertical bars indicate the frequencies uij — Vj 
of first order sideband resonances corresponding to 10 axial 
vibrational modes of 10 171 Yb + ions in an ion trap with the 
field gradient shown in Fig. |7| The frequencies are given 
relative to the transition frequency wi of the first ion, and the 
bars are labelled with j = 1, . . . , 10. The height of the bars 
indicates the transition probability associated with a given 
sideband relative to the first sideband (COM mode) of ion 1 
for a Raman transition. The relative transition probability is 
proportional to \Vj%\/\vli \ « Sfu~ (1/2) /Slu^ (1/2) with j = a 



resonances all lie within a frequency interval of about 
0.015 x ^i/2tt = 15kHz. 

The height of the bars in Fig. [S] indicates the strength 
of the coupling between the driving radiation and the 
respective sideband transition relative to the COM side- 
band of ion number 1 . The relevant coupling parameter is 
the LDP. For optical transitions \rf^ l/l^ffl ~ l^l/Wl = 
Sfv~ {1/2) /S{v~ (1/2] with j = a. The ratio of these pa- 
rameters for the highest vibrational mode [a = 10) is 
about 3 orders of magnitude smaller than for mode 1, 
since ion 10 is only slightly displaced from its equilibrium 
position when mode 10 is excited (compare the discussion 
in section ITlIBIl . 

We will now investigate how well simultaneous cooling 
can be done with the sideband resonances not perfectly 
superposed. 

B. Simultaneous Raman cooling with non-ideal 
gradient 

In Fig. |5Jl) the steady state vibrational excitation, 
(n'°')/ of a string of 10 ions is displayed as a function 
of the detuning of the Raman beams relative to the res- 
onance frequency u>i of ion 1. The Rabi frequencies and 
detuning, too, are the same as have been used to generate 
Fig. 0] However, the field gradient that shifts the ions' 
resonances is not assumed ideal as in Fig. 0] instead the 
one generated by three single windings as described above 
(Fig. 0) nas been used. Despite the imperfect superpo- 
sition of the cooling resonances, low temperatures of all 
modes close to their ground state can be achieved as can 
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FIG. 9: Raman sideband cooling a chain of 10 ions in the 
presence of a spatially inhomogeneous magnetic field. The 
condition cji — v\ — ... = ujn — vn is approximately ful- 
filled (see Fig. Same parameters as in Fig. 21 (i- e -j 
fii2 = 100 x 27rkHz, Q i = 5x 27rkHz, A i2 = -10 X 2ttMHz). 
a) Steady state vibrational excitation (n' a '}/ as a function 
of the detuning Aoi around Aoi = —v\. b) (rv- a ') f for each 
mode at that detuning where the sum of the mean vibrational 
quantum numbers of all ten modes is minimal. 



be seen in Fig. Eb). Here, the value of (n^ a ')f for each 
mode has been plotted at that detuning A i = — 1.008fi 
where the sum of all excitations is minimal. 

Fig. HUk displays the excitation of each mode over 
a wide range of the detuning such that all first order 
red sideband resonances are visible. Here, the Rabi fre- 
quency tt±2 of the repump laser has been increased to 
1 x 27rMHz as compared to 100 x 27rkHz in the previous 
figures. This results i) in higher final temperatures, and, 
ii) in broader resonances as is evident in Fig. 110b and 
thus makes cooling less susceptible to errors in the rela- 
tive detuning between laser light and ionic resonances. 

A higher steady state vibrational excitation due to the 
larger intensity of the repump laser is evident in Fig. 110b 
where (n^) f for each mode is plotted with the same 
parameters as in Fig. [UJa, however with H12 = 1 x 27rMHz. 

Errors and fluctuations in the relative detuning of the 
Raman laser beams driving the sideband transition are 
expected to be small and not to affect the efficiency of si- 
multaneous cooling, if the two light fields inducing the 
stimulated Raman process are derived from the same 




FIG. 10: Raman sideband cooling a chain of 10 ions in the 
presence of a spatially inhomogeneous magnetic field. The 
condition u>i — v\ = . . . = lum — um is approximately fulfilled 
(see Fig.|HJl. a) {rv-°^) / as a function of Aoi displaying all first 
order red sidebands with all parameters unchanged compared 
to Fig. |§] except O12 = 1 X 2ttMHz. b) {n {a) ) f as a func- 
tion of Aoi in the vicinity of the common sideband resonance 
around Aoi = —v\. The resonances are power broadened 
and ac-Stark shifted compared to the ones in Fig. Hi), c) 
Steady state population (n^ a ')f for each mode at that detun- 
ing where the sum of the mean vibrational quantum numbers 
of all ten modes is minimal (the analogue to Fig. [5Jj except 
that f2i2 = 1 X 27rMHz.) d) (n (o!) )/ for each mode at that 
detuning where one of the modes (here mode 10) reaches the 
absolute minimum. Ontima! cooline 1 of all vibrational modes 
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laser source using, for example, acousto-optic or electro- 
optic modulators. This is feasible by translating into 
the optical domain the microwave or radio frequency 
that characterizes the splitting of states |0) and |1). Mi- 
crowave or rf signals can be controlled with high precision 
and display low enough drift to ensure efficient cooling. If 
a large enough intensity of the repump laser is employed, 
then the steady state vibrational excitation varies slowly 
as a function of Aoi as is visible in Fig. 110b . Thus, 
the requirements regarding both the precision of adjust- 
ment and the drift of the source generating the Raman 
difference frequency are further relaxed. 

It should be noted that efficient cooling does not only 
occur around the resonance A i = —V\ but also at other 
values of A 01 as can be seen in Fig. HOfa . As an example, 
Fig. HBl shows (nW)/ of all modes at that detuning, 
Aoi = —3.91^i where mode 10 reaches its absolute min- 
imum. At this resonance the red sideband of the 5th 
ion corresponding to the 10th mode is driven by the Ra- 
man beams (compare Fig. Note that all other vibra- 
tional modes are also cooled at the same time. Therefore, 
an efficient procedure for cooling all vibrational modes 
close to their ground state would be to first tune the 
Raman beams such that mode 10 is optimally cooled 
(i.e., Aqi = —3.91^i, Fig. HUH), and subsequently set 
Aoi = — 1.022i/i (Fig. II Ob i in order to simultaneously 
cool modes 1 through 9. This approach is discussed in 
more detail in section ITlI B II 

Initial cooling of vibrational modes often is prerequi- 
site for subsequent coherent manipulation of internal and 
motional degrees of freedom of an ion chain, for example, 
quantum logic operations. When implementing quantum 
logic operations it may not be advantageous to address 
all motional sidebands with a single frequency as is done 
here for simultaneous sideband cooling. The magnetic 
field gradient that superposes the sideband resonances 
for cooling may then be turned off adiabatically after ini- 
tial Raman cooling of all vibrational modes. This should 
be done fast enough not to allow for appreciable heating 
of the ion string, for example, by patch fields, and slow 
enough not to excite vibrational modes in the process. 
A lower limit for the time it takes to ramp up the gra- 
dient seems to be 2ir/v\. Thus, with v\ = 1 X 27rMHz 
the additional time needed to change the field gradient is 
negligible compared to the time needed to sideband cool 
the ion string which for typical parameters takes between 
a few hundred /is and a few ms (compare Fig. |5j). 

If microwave radiation is used to coherently manip- 
ulate internal and motional degrees of freedom and for 
quantum logic operations, it is useful not to turn off the 
magnetic field gradient, but instead to ramp up the field 
gradient to a value where all coincidences between inter- 
nal and motional resonances are removed 0,^3- Also, a 
larger field gradient for quantum logic operations is de- 
sirable in this case to have stronger coupling between in- 
ternal and external states (i.e., a larger effective LDP rf^ 
[eq. I20|L The considerations in the previous paragraph 
regarding the time scale of change of the field gradient 



apply here, too. 

For the cooling scheme introduced here to work, the 
field gradient has to vary in the axial direction and it 
remains to be shown in what follows that this variation is 
compatible with the neglect of higher-order terms in the 
local displacement qj of ion j and in d z B in the derivation 
of the effective Lamb-Dicke parameter induced by the 
magnetic field gradient |l8|, LL9j . Neglecting higher order 
terms is justified as long as \q 2 jd 2 z B\ <§; \qjd z B\. Usingl^Bl 

and qj rs Az = ^/fij2mv\ this condition can be written 
as 



(Vj - Uj_i)(C: 



0) 



1 



Considering the region where the second derivative of the 
magentic field is maximal (J — 9) and inserting num- 
bers into relation E| gives for ten ions 0.24 <C 2.9 x 
10 3 (to/i^i)6. The right-hand side of this inequality is 
dimensionless if m is inserted in a.m.u., and yields « 500 
for 171 Yb + ions and V\ = 1 x 27rMHz. Hence, relation 1771 
is fulfilled. This can be understood by considering that 
the typical distance, Co over which the gradient has to 
vary is much larger than the range of motion, qj ss Az of 
an individual ion, and the approximation of a linear field 
gradient is a good one. 



VI. CONCLUSIONS 

We have proposed a scheme for cooling the vibrational 
motion of ions in a linear trap configuration. Axial vi- 
brational modes are simultaneously cooled close to their 
ground state by superimposing the red motional side- 
bands in the absorption spectrum of different ions such 
that the red sideband of each mode is excited when driv- 
ing an internal transition of the ions with monochromatic 
radiation. This spectral property is achieved by apply- 
ing a magnetic field gradient along the trap axis shifting 
individually the internal ionic resonances by a desired 
amount. 

Exemplary results of numerical simulations for the case 
of an ion chain consisting of N — 10 ions are presented 
and extensively discussed. Detailed simulations have also 
been carried out with 1 < N < 10 and for some values 
N > 10. They lead to the same qualitative conclusion as 
presented for the case of N = 10. 

Numerical studies show that simultaneously Raman 
cooling all axial modes is effective for realistic sets of pa- 
rameters. These studies also reveal that using microwave 
radiation to drive the sideband transition is not as effi- 
cient, due to the relatively small mechanical effect associ- 
ated with the excitation of this transition. The mechan- 
ical effect could be enhanced by applying a larger mag- 
netic field gradients. For simultaneous sideband cooling 
of all vibrational modes, however, the gradient is fixed by 
the requirement of superimposing the sidebands. Usual 
sideband cooling on a hyperfine or Zeeman transition us- 
ing microwave radiation becomes possible, if a larger field 
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gradient is used. This may be particularly useful, if a 
single ion (or an atom in an optical dipole trap) is to be 
sideband cooled using a microwave transition, or, if the 
vibrational modes of an ion chain are to be cooled sequen- 
tially using microwave radiation. Moreover, these tech- 
niques could also be implemented for sympathetic cooling 
of an ion chain, whereby some modes are simultaneously 
cooled by addres sing ions of other species embedded in 
the chain [S3, H3, E3| . 

In conclusion, we have shown that simultaneous side- 
band cooling with optical radiation can be efficiently im- 
plemented taking into account experimental conditions 
and even for a simple arrangement of magnetic field gen- 
erating elements. 
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VIII. APPENDIX 

In this appendix we introduce the hamiltonian and 
master equation describing the dynamics discussed in 
Sec. IIVBI We consider a chain consisting of N identical 
ions aligned along the z-axis, and in presence of a mag- 
netic field B(z). The internal electronic states of each ion 
which are relevant for the dynamics are the stable states 
|0) and |1) and the excited state |2). The transitions 
|0) — > |1), |1) — > 1 2) are respectively a magnetic and an 
optical dipole transition. We assume that the magnetic 
moments of |0) and |2) vanish, while |1) has magnetic 
moment fi. Thus its energy with respect to |0) is shifted 
proportionally to the field, ftiuo(z) cx |i?(z)|. The Hamil- 
tonian describing the internal degrees of freedom has the 
form: 

fl^ = nX)(wo(«i)|l>J<l|+wa|2>i<2|) (28) 



where the index j labels the ions along the chain. The 
collective excitations of the chain are described by the 
eigenmodes at frequency fx, . . . , vn, which are indepen- 
dent of the internal states. Denoting with Q a ,P a the 
normal coordinates and conjugate momenta of the oscil- 
lator at frequency v a , the Hamiltonian for the external 
degrees of freedom then has the form 



An 



N 



1 



2m ^ 



(29) 



+ - 



in 



N 



a = l 



2mv't 



E 



dzj 



\l)(l\S? 



where we have neglected the higher spatial derivatives of 
the magnetic field. 

Thus, the coupling of the excited state |1) to a spatially 
varying magnetic field shifts the center of the oscillators 
for the ions in the electronic excited state. 
The states |0) and |1) are coupled by radiation according 
to the Hamiltonian 



Woi 



E 

3 



1) 



(30) 



where the coupling can be generated either by microwave 
radiation driving the magnetic dipole, or by a pair of 
Raman lasers. In the first case, fioi is the Rabi frequency, 
woi the frequency of radiation and fcoi the corresponding 
wave vector. In the case of coupling by Raman lasers, 
VLqi is the effective Rabi frequency, u)q% the frequency and 
km the resulting wave vector describing the two-photon 
process. 

The transition |1) — > |2) is driven below saturation by 
a laser at Rabi frequency £li2, frequency W12, and wave 
vector k. The interaction term reads: 



Wv 



E 



|2) i (l|e _i ( u ' 12 *- fcz ^)e ife «J + H.c. 



(. 31 ) 

where qj is the displacement of the ion j from the classical 
equilibrium position Zj. 

The master equation for the density matrix p, describing 
the internal and external degrees of freedom of the ions, 
reads: 



Cp 



(32) 



where H = H + fl moc + Wqi + Wu, and Cp is the Li- 
ouvillian describing the spontaneous emission processes, 
i.e. the decay from the state |2) into the states |0) and 
|1) at rates T20, I^i, respectively, where T = T20 + T21 
is the total decay rate. The Liouville operator for the 
spontaneous decay is 



£p = -ir^[|2) i (2|p + p|2) i (2| 



(33) 



+r 20 / duAf(u)e iku *i 10), <2|p|2)j (0\ e - iku ^ 



+r 21 / duN{u)e lkuqj (2|p|2) j (1 |e~ ifetl<?J 



with N(u) the dipole pattern for spontaneous emission 
and u the projection of the direction of photon emission 
onto the trap axis. In order to study the dynamics, it 
is convenient to move to the inertial frames rotating at 
the field frequencies. Moreover, we apply the unitary 
transformation [19| 



U = exp 



2mv' 



E 



duj . 



Pa 



(34) 
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We denote with p the density matrix in the new reference 
frame. The master equation now reads 



(35) 



where H = Hq + H mcc + Wqi + W12, and the individual 
terms have the form: 

H = hJ2 [*(*i)|0>i<0| + A(z J )|2)(2|] (36) 



with S(zj) = luqi —uo(zj) and A(zj-) = u>% — u>o(zj) — uj\2- 
The mechanical energy is 



and rjj^ is an effective Lamb-Dicke parameter defined 
in Eq. (J2UJ. The dynamics of the individual modes are 
effectively decoupled from the others in the La mb-D icke 
regime, which holds when the condition (n a ) <C 1 

is fulfilled. 

The cooling rate and steady state occupation are eval- 
uated for each mode by using the procedure outlined 
in |30|. The rates A±, which determine the cooling 
rate and the steady state mean occupation according to 
Eqs. (0 and JHJl, are given by 



A% = 2Re{S a (Tv a )+D a } 



(42) 



Hmec — 



a—1 a—1 

^hv a (ala a + - J (37) 



where a^ a , a a are the creation and annihilation operator, 
respectively of a quantum of energy hv a . The interaction 
term between the states |0) and |1) now reads 



w 01 = J2 



Ml 



01 



|l)j(U|c 



\{k iZj—il>) 



ik 1 



k?P a 



where 



£a _ duJ Qj 1 

3 dzj 2mv^ 3 



H.c 

(38) 

(39) 



Thus, the excitation between two states, where the 
mechanical potential in one is shifted with respect to 
the other, corresponds to an effective recoil, here de- 
scribed by the effective Lamb-Dicke parameter ffj = 
k" ^jKmv a /2. Finally, the optical pumping between the 
states |1) and |2) is given by 

HQi2 



W l2 = ■ 



|2)j(l|e iXj + ^ Q * aQ) +ILc 



(40) 

where Xj is a constant phase, that depends on the posi- 
tion according to 



Xj 



k Zj -iP-hk^2kfSf/2 



(41) 



where Sa^fVa) is the fluctuation spectrum of the dipole 
force Fj a , 



S a {±v a ) = —— V / dte^TriFj^Fj^pst} 

IVl V a Jq 

(43) 

and D a is the diffusion coefficient due to spontaneous 
emission, 



D a = 0]T (|Wa| 2 r 20 + l?7i2 J j 2 r 2 i) (2|p st |2) (44) 

3 



where <f> = J duM(u)u 2 and here <j) — 2/5. Here, pst 
is the stationary solution of Eq. Q35J) at zero order in 
the Lamb-Dicke parameter, and the dipole force Fj, a is 
defined as 



Fj,a =i%i^|l)<0|+Wa^|2)(l|+H.c. (45) 



The fluctuation spectrum and the diffusion are found by 
evaluating numerically the steady state density matrix. 
The two-time correlation function in (|43|l is found by 
applying the quantum regression theorem according to 
the master equation (|35|l at zero order in the Lamb-Dicke 
parameter. 
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